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Using Monte Carlo simulations we have studied the transition from an "active" steady state to an ab- 
sorbing "inactive" state for two versions of the branching annihilating random walks with parity conservation 
on a square lattice. In the first model the randomly walking particles annihilate when they meet and the 
branching process creates two additional particles; in the second case we distinguish particles and antiparti- 
cles created and annihilated in pairs. Quite distinct critical behavior is found in the two cases, raising the 
question of what determines universality in this kind of systems. 
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Branching annihilating random walks (BARWs) have 
been extensively studied in recent years because they are 
a prototype for a variety of reaction-diffusion- like systems 
(for recent reviews see the work by Cardy and Tauber 
jl| and by Marro and Dickman Q for more general as- 
pects). For different models, the random walkers can 
represent either domain walls (kinks) or active sites on 
a lattice. In general, the corresponding critical behav- 
ior belongs to the directed percolation (DP) universality 
class. According to the "DP conjecture" ||, most of 
the one-component models with a single absorbing state 
belong to the DP universality class. Exceptions can ap- 
pear when additional symmetries or conservation laws 
are introduced. The best known one-dimensional (ID) 
exceptions are models in which, either the parity of the 
number of particles is conserved during the elementary 
processes ji] - || , or there is an equivalence between two 
absorbing states |l^,0. Henceforth, we will concentrate 
on this type of 2D BARWs whose behavior differs from 
the 2D DP universality class jjJJ . 

A field theoretical analysis of such systems was recently 
reported @,||. In ID, the numerical and theoretical ap- 
proaches are in satisfactory agreement: a new universal- 
ity class appears when parity is preserved. Agreement is 
also found for D > 2 in which case the mean-field results 
are valid. The nonexistence of the absorbing state (ex- 
tinction) for finite values of the branching rate predicted 
by mean-field approximation for D > 2 was also con- 
firmed by early Monte Carlo (MC) simulations |l4| in 2D 
and 3D. However, according to the field theoretical analy- 
sis, logarithmic corrections are expected for the marginal 
dimensionality D = 2 for models with parity conserva- 
tion and this was not born out by previous work. In this 
communication we report results of extensive MC simula- 
tions in 2D lattices showing evidence of such corrections. 
Inspired by the kink/antikink interpretation of creation 
and annihilation, we have also considered a modification 
of the model where two types of particles are present in 



the creation and annihilation processes and indeed found 
that a distinct behavior appears in this case. 

In the first model there is only one species of particles 
walking randomly on a square lattice. In addition to the 
single particle diffusion, the time evolution is governed by 
creation and annihilation of particle pairs as follows. A 
randomly chosen particle creates two additional particles 
with a probability p which are located on two randomly 
chosen neighboring empty sites fHI ; otherwise the chosen 
particle jumps to one of its nearest neighbor positions. 
For both elementary processes, if the destination site is 
occupied then the two particles annihilate. In order to 
study extinction, we consider only even initial numbers 
of particles. 

In the second model we distinguish particles and an- 
tiparticles, the creation and annihilation processes in- 
volving a particle/ antiparticle pair with an evolution rule 
similar to the above one. During the sequential updat- 
ing we neglect all the elementary processes which would 
result in two particles (or antiparticles) on the same site. 
The numbers of particles and antiparticles are chosen to 
be equal. As a result we have a unique absorbing state 
(no particles) independent of time. This model can be 
considered as a 2D generalization of the parity conserving 
model introduced by Menyhard J5| where the ID ferro- 
magnetic domain walls are represented by particles and 
antiparticles. It is worth mentioning that in the former 
model the particles and antiparticles are alternately po- 
sitioned along the chain and this feature is maintained 
by the elementary processes, therefore their distinction 
is not relevant. In the present 2D model, however, two 
particles (antiparticles) can occupy neighboring positions 
and they can avoid each other. 

The MC simulations are performed on an L x L square 
lattice with periodic boundary conditions for different 
values of the branching rate p. In order to have suffi- 
ciently accurate results, the system size (L) is increased 
up to L = 2000 for small concentration of particles. The 
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simulations are started from a randomly half-filled lattice 
and during the evolution we record the concentration of 
particles (other initial conditions were also tested). Time 
is measured in Monte Carlo steps (MCS) within which 
each particle has an opportunity to jump (probability 
1 — p) or branch (probability p) . 

For both models, in the absence of branching (p — 
0) , the number of particles decreases monotonously and 
eventually vanishes. For p > 0, however, the system 
remains active with a fluctuating number of particles n 
if the lattice size is sufficiently large. In the stationary 
state the average concentration (c = (n) / L 2 ) of particles 
vanishes continuously when p tends to 0. In other words, 
decreasing the branching rate, both systems undergo a 
transition with a critical point p c = 0. 

First we have investigated the decrease of concentra- 
tion at the critical point (p c — 0). For this purpose 
the time-dependent average concentration and its fluc- 
tuation x = L 2 {{c — n/L 2 ) 2 ) are determined at discrete 
time steps (equidistant in the logarithmic scale) by aver- 
aging over 500 runs. 
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FIG. 1. Concentration (upper curve) and its fluctuation 
(lower curve) multiplied by the time vs. the logarithm of 
time at zero branching rate for the first model. 

For the first model the field theoretical investigation 
suggests c oc hit/t Q. To check this prediction we have 
plotted tc vs. t in a log-lin plot. Our MC results (see 
Fig. |l|) obtained for L = 2000 indicate clearly that the 
time-dependent concentration can be well described as 



for sufficiently long times (t > 100 MCS). This function 
fits the MC data if A = 0.2238 and B = 0.8979. The 
leading term of this asymptotic behavior agrees with the 
above prediction given by Lee |]l6| and Cardy and Tauber 

The fluctuation of concentration decreases proportion- 
ally with the average value of concentration as demon- 
strated in Fig. |]. Neglecting the "noisy decoration" due 



to the statistical error, the ratio of the time-dependent 
fluctuation and concentration can be well approximated 
as x{t)/ c (t) — 0.67(2) on the time range indicated in 
Fig. W. This ratio agrees very well with the theoretical 
prediction (2/3) obtained by Lee jl6| using renormaliza- 
tion group technique. 

For the second model the decrease of concentration 
follows a different behavior at the critical point. Indeed 
when p — we get the diffusion limited surface reac- 
tion A+B — > case, already studied by several authors 
jl7]-[n|. The concentrations decrease as t~ d / 4 if the di- 
mension d is lower than 4 (the upper critical dimension 
for this system). Compared to the former case, the re- 
sults depend more strikingly on size, as shown in a log-log 
plot (Fig. H). The sharp decrease in c is a consequence of 
the extinction whose probability is higher for smaller sys- 
tems. Our simulation confirms the mentioned power law 
behavior in the limit L — > oo as indicated in Fig. 0. The 
extent of this behavior for the largest system is illustrated 
in the figure by the dashed line. The fluctuations are al- 
most constant during this scaling regime (x(t) — 0.009(1) 
for L = 2000 and 600 < t < 60000). 




FIG. 2. Time dependence of concentration of particles (and 
antiparticles) in the second model at zero branching rate for 
different system sizes (L = 100, 200, 300, 500, 1000 and 2000 
from left to right). The dashed line (slope —0.5) indicates the 
theoretical power law behavior. 

Systematic and extensive MC simulations have been 
performed to study the average concentration and its 
fluctuation in the steady state (reached after some ther- 
malization) for finite branching rates. For the smallest p 
values both the thermalization and sampling times were 
longer than 10 5 (10 6 ) MCS for the first (second) model. 
These simulations were repeated 20 times to suppress the 
undesired effects of long time fluctuations. 

The results for the first model are summarized in 
Fig. ||. At first glance, the MC data for the concen- 
tration (diamonds in the log-log plot) indicates a power 
law behavior, namely c(p) oc p@ with (3 = 1.276, similar 
to what was found by Takayasu and Tretyakov |H| . The 
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careful reader can, however, observe a definite deviation 
from this behavior (positive curvature) whose magnitude 
exceeds our statistical error. 




FIG. 3. Log-log plot of the particle concentration (dia- 
monds) and its fluctuation (squares) as a function of branch- 
ing rate for the first model. The solid line indicates the fitted 
curve described in the text. 

Taking the logarithmic corrections into account Cardy 
and Tauber have suggested |D| that the leading term 
of the c(j>) function is proportional to p/hi 2 (p). This 
function does not fit adequatly the present MC data; 
however, an excellent fitting is found if we use c(p) = 
p/[A + Sln(p) + Cln 2 (p)]. This function is represented 
by a solid line in Fig. | for A = 0.4098, B = -0.5825 
and C — 0.09601. Notice that this formula confirms the 
theoretical prediction in the limit p — > 0. 




FIG. 4. Particle and antiparticle concentration (diamonds) 
and its fluctuation (squares) vs. branching rate for the second 
model. The fitted asymptotic power law behavior is indicated 
by the straight line (slope 1.11). 

Figure || indicates clearly that the concentration fluc- 
tuation is proportional to the concentration itself in the 



p region we have studied. These quantities satisfy the 
relation x(p)/ c {p) — 2.5(2) within the statistical error. 

Similar investigations have been performed for the sec- 
ond model. Figure [| demonstrates clearly that the MC 
data tends towards a power law for both the concentra- 
tion (c) and its fluctuation (x) a ^ small values of branch- 
ing p. Fitting the function c(p) = Ap@ to our numerical 
data we have obtained that (3 = 1.11(1). The ratio of the 
fluctuation to the concentration is smaller than found for 
the first model, namely x(p)/ c (p) = 1-10(5). 

In summary, we have studied and compared two simple 
models of BARWs on a square lattice with parity conser- 
vation. In contrast with the first model - where there is 
only one type of particles - the second model has parti- 
cles and antiparticles annihilating (only) each other when 
they meet. This distinction has caused significant differ- 
ences between their behavior at the critical point (no 
branching) as well as in the stationary states for finite 
branching rate. For the first model, our MC simulations 
have justified the appearance of the logarithmic correc- 
tions predicted theoretically by Lee and by Cardy and 
Tauber. On the contrary, we have observed power law 
behavior in the second model. Surprisingly, the fluctua- 
tions decrease with p - more precisely, the p dependence 
of the fluctuation is found to be proportional to the con- 
centration for both models. The significant differences 
between the behavior of the present models imply the 
possibility to find other two-dimensional systems whose 
critical behavior does not belong to the directed percola- 
tion universality class. The belief that in these systems 
parity conservation is sufficient to determine the univer- 
sality class is probably also to be questioned. 
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